Effect of electron-doping on spin excitations of underdoped BaFei 96Nio.o4As2 
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We use neutron scattering to study magnetic order and spin excitations in BaFei.96Nio.o4As2. On 
cooling, the system first changes the lattice symmetry from tetragonal to orthoromhbic near ~97 
K, and then orders antiferromagnetically at Tjv = 91 K before developing weak superconductivity 
below ~15 K. Although superconductivity appears to co-exist with static antiferromagnetic order 
from transport and neutron diffraction measurement, inelastic neutron scattering experiments reveal 
that magnetic excitations do not respond to superconductivity. Instead, the effect of electron- 
doping is to reduce the c-axis exchange coupling in BaFe2As2 and induce quasi two-dimensional 
spin excitations. These results suggest that transition from three-dimensional spin waves to two- 
dimensional spin excitations by electron-doping is important for the separated structural/magnetic 
phase transitions and high-temperature superconductivity in iron arsenides. 
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Antiferromagnetism is relevant to high temperature 
(high- Tc) superconductivity in copper oxides and iron ar- 
senides because superconductivity arises from electron- 
or hole-doping of their static antiferromagnetic (AF) or- 
dered parent compounds P, i, i, i, i, 0, 0, S H- In 
the case of cuprates, spin waves of the parent materials 
can be very well described by a local moment Heisen- 
berg Hamiltonian and spin excitations in optimally doped 
superconductors are dominated by a neutron spin reso- 
nance centered at the AF ordering wavevector [J, 0, [3] . 
For undoped iron arsenides such as AFegAsg {A =Ba, 
Sr, Ca) with a spin structure of Fig. la [lO|, spin waves 
consist of a large anisotropy gap at the AF zone cen- 
ter [A(l, 0, 1) < 9.8 me VI and excitations extend up to 
-200 meV [uljll (Tl, [ij, [ll| . For optimally doped su- 
perconductors [Sl, S III J the gapped spin wave excitations 
were replaced by a gapless continuum of scattering in 
the normal state and a neutron spin resonance below Tc 
Tdl . 17, 1^ Si- Since spin fluctuations may play a cru- 



the superconductivity 28|, |29|] . Although recent neutron 
scattering experiments confirmed that the upper transi- 
tion is structural and the AF order occurs at lower tem- 
perature [sol . 31 . [s^l, its microscopic origin is still un- 



cial role in the electron pairing and superconductivity of 
iron arsenides 2^, 21, 22, 23 , [ij, [23, [26 , ^] , it is impera- 
tive to determine how the spin dynamics of the undoped 
AF parent compounds evolve as they are tuned toward 
optimally doped superconductivity by electron or hole 
doping. 

In the undoped state, BaFe2As2 exhibits simultaneous 
structural and magnetic phase transitions below Tg = 
T/v = 143 K, changing the crystal lattice symmetry from 
the high-temperature tetragonal to low-temperature or- 
thorhombic phase [13] . Upon Co-doping to induce elec- 
trons onto the FeAs plane, the combined AF and struc- 
tural phase transitions were split into two distinct tran- 
sitions and the electronic phase diagram in the lower Co- 
doping region displays coexisting static AF order with 



known. More importantly, it is unclear what happens to 
the spin waves of BaFe2As2 when electrons are doped into 
these materials. While Pratt et al. [s^ reported gapless 
normal state spin excitations for BaFei.906Coo.o94As2, 
measurements on BaFei .AzCon nsAsa suggest gapped nor- 
mal state excitations [32|. On cooling below Tc, both 
materials reveal a reduction in the static ordered AF mo- 
ment and the appearance of a spin resonance |30l . i32,] . 

To compare with the results obtained on Co-doped 
BaFe2As2 [H, [i^, [13, 31, [13|, we carried out neutron 
scattering experiments on Ni-doped BaFei.96Nio.o4As2 
(Tc « 15 K, Figs. lc,ld) [TJ- In contrast to the results on 
Co-doped materials 30|, [S^ , we find that the static AF 
order and spin excitations in BaFei. 96Nio.o4As2 do not re- 
spond to the occurence of superconductivity. Instead, the 
effect of electron-doping is to significantly reduce the c- 
axis exchange coupling and change the three-dimensional 
(3D) spin waves of BaFe2As2 into quasi two-dimensional 
(2D). These results suggest that the separated struc- 
tural/magnetic phase transition and the appearance of 
bulk superconductivity upon doping may be associated 
with the diminishing spin anisotropy gap and the 3D to 
2D transition of the spin excitations. 

Using the self-flux method 7] , we grew a 1 gram sin- 
gle crystal of BaFei.g6Nio.o4As2 with an in-plane and out- 
of-plane mosaic of 1.74° and 2.20° full- width at half maxi- 
mum (FWHM), respectively. We defined the wave vector 
Q at (ga;, qy, g^) as {H,K,L) = {q^a/2TT,qyb/2TT,q^c/2TT) 
reciprocal lattice units (rlu) using the orthorhombic mag- 
netic unit cell (space group Fmmm), where a = 5.5 A, 



2 



ib) 



BaFe, Ni As. 



\ • Bulk Superconductivity 
V ■ Filamentary Superconductivity 



;=o.oi 



G_ 0.04 ■ 
0.02 ■ 


C) T - 97K 

=' — ^ 

: t 

• Tj,= 91K 


0.000 ■ 
^ -0.004 ■ 


0.00 ■ 


j 






100 2 

T(K) 





0.0035 ■ 
■§ 0.0030 ■ 


»)| 1 1 (2,0,0) 




5 0.0025 ■ 


T Cooling 

Ts~97K^ + 1 

Tn = 91K I 1 


Amplitude | 


T(K) 



/ 



• Zero Field Cooled (5 Oe} 
■ Field cooled (5 Ofl) 




3) • I T = 18KData 
A I T = 18K BKG 
■ I T = 3K Data 
▼ I T = 3K BKG 1 

13 



L * T - 100K BKG 




T T = 1D0K 
♦ T = 85.6K 

.'. T = 60K 



10 2 



B 8 10 



E (meV) 
H (rlu) 




FIG. 1: (color online), (a) Diagram of the parent com- 
pound BaFe2As2 with Fe spin ordering and magnetic ex- 
change couplings depicted. We use the same unit cell for 
BaFei.g6Nio.o4As2. (b) Electronic phase diagram from Ref. 
01 • (c) Temperature dependence of the resistance showing 
anomalies at Ts, Tm, and T^. (d) Temperature dependence 
of the Meissner and shielding signals on a small crystal (field 
cooled Atix = —0.001 at 4.5 K). (e) Temperature dependence 
of the structural distortion of the lattice as determined by 
tracking the width of the (2, 0, 0) nuclear Bragg peak using 
A/2 scattering without Be filter, (f) Magnetic order param- 
eter determined by Q-scans around (1,0,1) magnetic Bragg 
peak above background. The solid line shows order parame- 
ter fit using oc (1 - T/TNf^ with Tn = 91.3 ± 0.7 K and 
/3 = 0.3 ±0.02. 



b ~ 5.4 A, and c = 12.77 A. We performed our neutron 
scattering experiment on the PANDA cold triple-axis 
spectrometer at the Forschungsneutronenquelle Heinz 
Maier-Leibnitz (FRM II), TU Munchen, Germany as 
described earUer [l^. Our sample was aligned in the 
[H, 0, L] zone inside a closed cycle refrigerator. 

Figures Ic and Id show the resistivity and susceptibil- 
ity data. The resistivity shows clear anomalies near 97 
K and 91 K before superconductivity sets in below '^15 
K (Fig. Ic). Although the presence of superconductiv- 
ity below Tc « 15 K is confirmed in the susceptibility 
measurement (Fig. Id), the weak Meissner effect sug- 
gests superconducting volume fraction of less than 0.2%. 



FIG. 2: (color online), (a) Energy scans at Q = (1,0,1) 
and Q = (1.2,0,1) above and below Tc. (b) x'iQy^) s-t 
Q = (1,0,1). (c) Energy scans at higher temperatures and, 
(d) the corresponding x" iQy {^) Q-scans along the [H, 0, 1] 
direction at 4 meV and different temperatures. Fourier trans- 
forms of the Gaussian peaks give minimum correlation lengths 
of C = 57±4 A(f) Estimated x"(Q,w) at 4 meV. (g) x"(Q,'^) 
at 7 meV with minimum correlation length of ^ = 54±6 A. (h) 
Low temperature Q-scans along the [1, 0, L] direction (c-axis) 
at 4 meV and 7 meV. The minimum dynamic spin correla- 
tion lengths are ^ ~ 14 ± 5 and 21 ± 2 A for 4 and 7 meV, 
respectively. 



Similar to Co-doped BaFe2As2 3^, 31, 32 1, we find that 
the tetragonal to orthorhombic structural transition hap- 
pens at 97 K while the AF order occurs below = 91 
K (Figs, le and If). However, in contrast with Co- 
doped materials 3^, 32 1, superconductivity has no influ- 
ence on the static AF order of BaFei.96Nio.o4As2 (Fig. 
If). This is consistent with the fact that our sample has 
a lower electron-doping to the FeAs-plane than those of 
33,111. 

the undoped BaFe2As2, spin waves have an 
anisotropy gap about 8 meV at Q = (1, 0, 1) [A(l, 0, 1) = 



Refs. 
In 



8 meV] [1^ uM- For optimally Co and Ni doped ma- 



terials, spin excitations are gapless in the normal state 
[itI . [lij and superconductivity-induced spin gaps open 
below Tc [l^. Figure 2a shows the constant-Q scans at 
the Q = (1, 0, 1) (signal) and Q = (1.2, 0, 1) (background) 




FIG. 3: (color online), (a) Energy scans at Q — (1,0,0) and 
Q = (1.4, 0, 0) from 0.5 meV to 7 meV at 3.5 K and 18 K. (b) 
Background corrected x"{Qj^) showing clear evidence for a 
4 meV spin gap, larger than that at Q = (1,0,0). (c) Tem- 
perature dependence of the signal [Q — (1,0,0)] and back- 
ground [Q = (1.4,0,0)] scattering at various temperatures, 
(d) x'iQy'^) at 4 meV. (e) Q-scans along the [i/, 0,0] direc- 
tion at 4 meV and different temperatures. A peak centered at 
Q = (1,0,0) appears above 80 K. (f) Background corrected 
x"{Q,^)- (g) Temperature dependence of the Q-scans along 
the [H, 0, 0] direction at 6 meV. The scattering has a peak at 
3 K. (h) Temperature dependence of the x"{Qj^) a-t 6 meV. 



positions above and below Tc for BaFei.96Nio.o4As2. Fig- 
ure 2b pfots the imaginary part of the dynamic suscepti- 
bihty x"(Q,w) after correcting for background and Bose 
population factor. We find that x"{Qi^) has a 2 meV 
spin gap and is not affected by superconductivity. Fig- 
ures 2c and 2d reveal that the magnetic intensity increase 
with increasing temperature below Tjv is due mostly to 
the Bose population factor. These results are confirmed 
by Q-scans along the [H, 0, 1] direction at different tem- 
peratures (Figs. 2e-2g), which display well-defined peaks 
at Q — (1,0,1) that have similar widths to the undoped 
BaFe2As2 at 10 meV [l^. Figure 4h shows Q-scans 
along the c-axis [1,0, L] direction. In contrast to the in- 
plane spin excitations, the c-axis spin-spin correlations 
are much broader for BaFei.96Nio.o4As2, thus suggesting 
2D nature of the excitations. 



FIG. 4: (color online), (a) Temperature dependence of the 
1 meV scattering at the signal Q = (1,0,0) and background 
Q = (1.4,0,0) positions. The inset shows Q-scans along the 
[H, 0, 0] at 1 meV and different temperatures. The scattering 
shows no anomaly across but clearly peak at Tm- (b) 
Temperature dependence of the scattering at 4 meV and Q = 
(1,0, 1) again peaks at Tjv. 



Further evidences for 2D spin excitations in 
BaFei.96Nio.o4As2 are summarized in Fig. 3. As- 
suming spin excitations in BaFe2-a;Nia;As2 can be 
described by an effective Heisenberg Hamiltonian, the 
spin anisotropy gaps at Q = (1,0, 1) and Q = (1, 0, 0) are 

A{1,0,1) = 2S[iJia + 2J2 + Jc + Js?-{Jc + Jla + 2J2f]'/^ 

and A(1,0,0) = 2SU2Jia + 4J2 + J.)(2J, + J,)]i/2, 
respectively [ll|, [13, 14, 15 1. Here S is the magnetic 
spin (= 1); Jia, J2, Jc, Js are effective in-plane 
nearest-neighbor, next nearest-neighbor, c-axis, and 
magnetic single ion anisotropy couplings, respectively 
(Fig. la). For BaFe2As2, we estimate A(1,0, 1) = 7.8 
meV and A(1,0, 0) = 20.2 meV assuming Jiq = 36, 



12, 14, n 



Upon 



J2 = 18, Jc = 0.3, Js = 0.106 meV 
electron doping to form BaFei.96Nio.o4As2, these spin 
gap values have been reduced to A(1,0, 1) — 2 meV 
and A(1,0,0) = 4 meV (Figs. 2b and 3b). Since such 
electron-doping hardly changes the in-plane Q-scan 
widths compared to that of the undoped BaFe2As2 
(Figs. 2e-g, 3e,3g) 'l^, it should only slightly modify 
the in-plane exchange couplings. Assuming that Jia and 
J2 are unchanged in BaFei.96Nio.o4As2, the observed 
A(l, 0, 1) = 2 meV and A(l, 0, 0) = 4 meV would corre- 
spond to Jc — 0.01 meV and Js — 0.007 meV, suggesting 
a rapid suppression of c-axis exchange coupling and 
magnetic single ion anisotropy with electron doping. 
In Ref. [3^ , it was argued that spin anisotropy for 
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BaFei.g2Coo. 08^-82 is similar to that of the BaFe2As2, 
meaning that the reduction in spin gap at Q = (1,0, 1) 
arises mostly from reduced Jia and J2- Assuming the 
best fitted values of S{Jia + 2 J2) = 32 meV and SJ^ — 
0.34 meV HH, we expect A(1,0,1) = 5.5 meV and 
A(1,0,0) = 14.2 meV with SJs = 0.106 meV. These 
values are clearly different from the observation. Even if 
we assume all exchange couplings to reduce by 50% upon 
electron-doping with S'(Jia + 2J2) = 32 meV, SJc = 0.15 
meV, and SJs = 0.05 meV, we still find A(1,0, 1) = 3.8 
and A(l, 0, 0) = 10 meV. This suggests that the large re- 
duction in the A(l, 0, 0) gap values upon electron doping 
is due to the reduced Jc and three-dimensionality of the 
system. 

To determine the temperature dependence of 
A(1,0, 0), we show in Fig. 3c the observed scat- 
tering at the signal Q — (1,0,0) and background 
(1.4,0,0) positions at several temperatures. Figure 
3d plots the estimated x"(Q,a;). Comparing Fig. 3d 
with Fig. 3b, the 4 meV spin gap x"{Qt^) at 18 K 
vanishes upon warming to above 60 K. These results 
are confirmed by Q-scans at 4 meV along the [ff, 0,0] 
direction (Fig. 3e). While scans at 2 K and 18 K are 
featureless, the scattering at 86 K and 100 K shows 
clear peaks centered at Q — (1,0,0). For Q-scans at 
6 meV, the scattering shows well-defined peaks at all 
temperatures (Fig. 3g). Converting these data into 
x"{Q,oj) in Fig. 3h confirms the results of Fig. 3d. 

Finally, we show in Fig. 4a the temperature depen- 
dence of the 1 meV scattering at the Q = (1,0,0) (sig- 
nal) and Q = (1.4,0,0) (background) positions. While 
the background scattering only increases slightly with in- 
creasing temperature and shows no anomaly across Tjv, 
the scattering at Q = (1,0,0) clearly peaks at T^r. Q- 
scans along the [H, 0, 0] direction at 1 meV confirm these 
results (the inset of Fig. 4a). Temperature dependence 
of the scattering at 4 meV and Q ~ (1,0,1) show similar 
behavior (Fig. 4b). These results suggest that the disap- 
pearing A(l, 0, 1) and A(l, 0, 0) gaps near arise from 
critical scattering associated with the static AF order. 

To understand the separated structural and magnetic 
phase transitions for BaFei q6 Nin.n 4As2, we note that in 
an effective J1-J2-JC model [22,l23[, the separation of the 
lattice and magnetic transition temperatures is controlled 
by the value of Jc 2^- There is only one transition tem- 
perature when Jc is large. A finite separation between 
the two transition temperatures occurs when Jcj Ji is 
reduced to the order of 10~^. Our experimental result 
of Jc/ J2 ~ 0.5 X 10"'^ is consistent with this picture. To 
quantitatively estimate the reduced T/v due to the smaler 
Jc, we note that Tjv ~ J2/ln{J2/ Jc) HI]. Let J° be the 
magnetic exchange values for the parental compounds, 
we can write Tpf/T^ — a[ln{a) + ln{b) + ln{c)]/ln{c), 
where a = J2/J2,b ~ Jc/Jc,c = J^/Jc- Using the ex- 
perimental values of the exchange coupling parameters 
determined earlier, we obtain (J2/J2) = (1/a) ^ 0.87, 



which is self-consistent with our suggestion that upon 
doping, the coupling between layers Jc is dramatically 
reduced while the change of the in-plane magnetic ex- 
change coupling is small. These results provide a natural 
and consistent interpretation for our experimental obser- 
vations. 

In summary, we have shown that the most dramatic 
effect of electron-doping in BaFe2As2 is to transform the 
3D anisotropic spin waves into 2D spin excitations. These 
results suggest that reduced dimensionality in spin exci- 
tations of iron asenides is important for the separated 
structural/magnetic phase transition and the occurance 
of high- Tc superconductivity in these materials. 
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